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ABSTRACT 


The orbital debris environment model contained in this report is intended 
to be used by the spacecraft community for the design and operation of 
spacecraft in low Earth orbit. This environment, when combined with 
material -dependent impact tests and spacecraft failure analysis, is 
intended to be used to evaluate spacecraft vulnerability, reliability, 
and shielding requirements. The environment represents a compromise 
between existing data to measure the environment, modeling of this data 
to predict the future environment, the uncertainty in both measurements 
and modeling, and the need to describe the environment so that various 
options concerning spacecraft design and operations can be easily 
evaluated. 


INTRODUCTION 

The natural meteoroid environment has historically been a design 
consideration for spacecraft. Meteoroids are part of the interplanetary 
environment and sweep through Earth orbital space at an average speed of 
20 km/s . At any one instant, a total of 200 kg of meteoroid mass is 
within 2000 km of the Earth's surface. Most of this mass is concentrated 
in 0.1 mm meteoroids. 

Within this same 2000 km above the Earth's surface, however, is an 
estimated 3,000,000 kg of man-made orbiting objects. These objects are 
in mostly high inclination orbits and sweep past one another at an 
average speed of 10 km/s. Most of this mass is concentrated in 
approximately 3000 spent rocket stages, inactive payloads, and a few 
active payloads. A smaller amount of mass, approximately 40,000 kg, is 
in the remaining 4000 objects currently being tracked by U.S. Space 
Command radars. Most of these objects are the result of more than 90 on- 
orbit satellite fragmentations. Recent ground telescope measurements of 
orbiting debris combined with analysis of hypervelocity impact pits on 
the returned surfaces of the Solar Maximum satellite indicate a total 
mass of approximately 1000 kg for orbital debris sizes of 1 cm or 
smaller, and approximately 300 kg for orbital debris smaller than 1 mm. 
This distribution of mass and relative velocity is sufficient to cause 
the orbital debris environment to be more hazardous than the meteoroid 
environment to most spacecraft operating in Earth orbit below 2000 km 
altitude . 

Mathematical modeling of this distribution of orbital debris predicts 
that collisional fragmentation will cause the amount of mass in the 1 cm 
and smaller size range to grow at twice the rate as the accumulation of 
total mass in Earth orbit. Over the past 10 years, this accumulation has 
increased at an average rate of 5 percent per year, indicating that the 
small sizes should be expected to increase at 10 percent per year. 
Reasons that both of these rates could be either higher or lower, as well 
as other uncertainties in the current and projected environment, are 
discussed in the section "Uncertainty in Debris Flux". As new data 
become available, a new environment will be issued. 

The authors wish to acknowledge those responsible for original research 
and data analyses utilized by this work. Dr. Andrew Potter, Mr. John 
Stanley, Dr. Karl Henize, Dr. Faith Vilas and Mr. Eugene Stansbery (all 


1 


DATA SOURCES 


The following data sources were considered in the construction of this 
environmental model: 

1. Orbital element sets supplied by US Space Command (both the 
cataloged population and those objects awaiting cataloging) 
for the period between 1976 and 1988. 

2. Optical measurements by MIT in 1984 using the telescopes of 
their Experimental Test Site (ETS) in Socorro, NM. 

3. Measurements designed to determine orbital debris particle 
albedo using a ground-based IR telescope at AMOS/MOTIF, US 
Space Command radars, and both NASA and Space Command 
telescopes . 

4. Analysis of hypervelocity impacts on the surfaces returned 
by the Shuttle from the repaired Solar Maximum Mission 
satellite in 1984. 

5. Mathematical models which consider various traffic models 
and satellite fragmentation processes to predict the future 
accumulation of debris . 


KEY ASSUMPTIONS AND CONCLUSIONS CONCERNING DATA SOURCES 


The following assumptions and/or conclusions were made or reached 
concerning the above data sources : 

1. The flux resulting from the US Space Command orbital 
element sets is complete to a limiting size of 10 cm for 
objects detected below 1000 kra altitude. 

2. The MIT telescopes observed a flux which is 5 times the 
flux predicted by US Space Command orbital element data 
sets . 

3. The MIT telescopes were observing objects to a limiting 
size of 2 cm in diameter (16 th magnitude at an albedo of 
0 . 1 ). 

4. The surfaces of the Solar Maximum Mission satellite 
experienced an orbital debris flux which varies from 20% of 
the meteoroid flux for debris sizes larger than 0.05 cm to 
a factor of 1000 times the meteoroid flux for sizes larger 
than 1 /xm. 

5. The orbital debris flux between 0.05 cm and 2 cm can be 
obtained by a linear interpolation (on a log 10 F (flux) vs 
log 10 ^ (diameter) plot) of the Solar Maximum Mission 
satellite surface data and the MIT telescope data. 



6. For any given size of orbital debris, the variation of flux 
with altitude, solar activity, orbital inclination, and the 
velocity and direction distribution is the same as that 
predicted by the US Space Command orbital element set data. 

7. The accumulation of objects tracked by US Space Command, 
when averaged over an 11-year solar cycle, will increase at 
a rate of 5% per year. 

8. The accumulation of objects detected by the MIT telescopes 
and the Solar Maximum satellite surfaces, when averaged 
over an 11-year solar cycle, will increase at twice the 
rate of the tracked objects, or 10% per year. 


DESIGN STANDARD 

I. Recommended Flux for Orbital Debris 

The cumulative flux of orbital debris of size d and larger on 
spacecraft orbiting at altitude h, inclination i, in the year t, when 
the solar activity for the previous year is S , is given by the 
following equation: 

F(d,h,i,t,S) = k^(h,S)^(i)*[F 1 (d).g 1 (t) + F 2 (d) • g 2 ( t) ] (1) 


where 

F - flux in impacts per square meter of surface area per year 
k * 1 for a randomly tumbling surface; must be calculated for a 
directional surface 
d = orbital debris diameter in cm 
t - time expressed in years 
h - altitude in km (h < 2000 km) 

S — 13 month smoothed solar flux F^q . 7 expressed in 10 4 Jy; 

retarded by one year from t 
i = inclination in degrees 

and 

*(h,S) = <^ 1 (h,S)/(^ 1 (h,S) + 1) 

<^ 1 (h,S) - 10< h / 200 ‘ s / 1 * 0 - 1 - 3 ) 

Fi<d) = 1.05x10' 5 *d“ 2 • 5 
F 2 (d) = 7 . OxlO 1 0 • (d + 700) " 6 

p, the assumed annual growth rate of mass in orbit =0.05 
gl (t) - (1 + 2 -p) (t * 1905 > 


g 2 (t) - (1 + p) (t ' 1905 > 

The inclination- dependent function f is a ratio of the flux on a 
spacecraft in an orbit of inclination i to that flux incident on a 



spacecraft in the current population's average inclination of about 
60°. Values for ip are given in figure 1 and tabulated in Table I. 

An average 11 -year solar cycle has values of S which range from 70 at 
solar minimum to 150 at solar maximum. However, the current cycle, 
which peaks in the year 1990, is predicted to be above average, 
possibly exceeding 200. 

An example orbital debris flux is compared with the meteoroid flux 
from NASA SP8013 in figure 2 for h - 500 km, t - 1995, k = 1.0, i- 
30°, and S = 90.0. 

The flux is defined such that the average number of impacts N on a 
spacecraft surface area of A exposed to the environment for the 
interval to t f is given by the following equation: 

N - f F-A dt (2) 

where A is the surface area exposed to the flux F at time t. 

The value of k can theoretically range from 0 to 4 (a value of 4 can 
only be achieved when a surface normal vector is oriented in the 
direction of a monodirectional flux) , and depends on the orientation 
of A with respect to the Earth and the spacecraft velocity vector. If 
the surface is randomly oriented, then k - 1. If the surface is 
oriented, with respect to the Earth then section IV must be used to 
calculate a value for k. In general, if the surface area is facing in 
the negative velocity direction, k - 0. However, if this area is 
facing in the same direction as the spacecraft velocity vector, and 
the spacecraft orbital inclination is near polar (which causes more 
"head-on" collisions), then k will approach its maximum value of about 
3.5 for the current directional distribution. 


The probability of exactly n impacts occurring on a surface is found 
from Poisson statistics, or 


P 

n 



-N 

e 


(3) 


II. Uncertainty in Debris Flux 

Factors which contribute significantly to the uncertainty in the 
orbital debris environment are inadequate measurements, an uncertainty 
in the level of future space activities, and the statistical character 
of major debris sources. The environment has been adequately measured 
by ground radars for orbital debris sizes larger than 10 cm. A 
limited amount of data using ground telescopes has shown a 2 cm flux 
which is currently estimated to be known within a factor of 3. 
Orbital debris sizes smaller than .05 cm have only been measured at 
500 km; at this altitude and for these smaller sizes, the environment 
is know within a factor of 2. Interpolation was used to obtain the 
flux between 0.05 cm and 2 cm at 500 km, and would be justified if the 
amount of mass between these two sizes were about the same as the mass 


contributing to the two sizes, or about 100 kg to 1000 kg. 
Mathematical modeling of various types of satellite breakups in Earth 
orbit make such an assumption seem reasonable. However, other than 
"reasonableness", there is no data which would prevent the flux of any 
particle in the size range between 0.05 cm and 2 cm to be as high as 
the 0.05 cm flux, or as low as the 2 cm flux, that is, vary by as much 
as several orders of magnitude. 

An additional uncertainty from the measurements arises because there 
are no measurements of debris smaller than 2 cm at other than 500 km 
altitude. Mathematical modeling concludes that if the debris is in 
near circular orbits and the source of the debris is at higher 
altitudes , the ratio of the amount of small debris to large debris 
should decrease with decreasing altitude. This ratio is assumed 
constant in the design environment. Consequently, there would be a 
smaller flux of less than 2 cm debris at altitudes less than 500 km, 
and a larger flux at altitudes above 500 km than is predicted by this 
model. However, if the debris is in highly elliptical orbits, then 
the flux of small debris could be nearly independent of altitude. 
Consequently, the amount that the flux differs from the design 
environment could be as high as a factor of 10 (either higher or 
lower) for every 200 km away from the 500 km altitude, up to an 
altitude of about 700 km. The large number of breakups at altitudes 
between 700 km and 1000 km and at 1500 km, together with the extremely 
long orbital lifetimes of fragments in these regions, make any 
predictions very sensitive to the nature of each of these breakups. 
The US Space Command data gives fluxes at 800 km and 1000 km which are 
twice as high as predicted by the recommended flux model, as shown in 
figure 3. For most altitudes between 1000 km and 2000 km, the current 
flux from objects tracked by US Space Command is significantly lower 
than the design environment. However, the large number of breakups at 
1500 km could have scattered smaller fragments over this region; in 
addition, future traffic may increase the flux of larger objects. 

Predicting future activity in space is highly uncertain. Since 1966, 
the non-US launch rate has increased by a average of 10% per year; 
however, US launch rates have decreased at this same rate, leading to 
a constant world launch rate since 1966. This constant launch rate 
has lead to a decreasing percentage growth in the accumulation of 
objects being tracked by the US Space Command. Averaged over the last 
solar cycle, this accumulation has grown at an average rate of 5% per 
year. A continued constant launch rate would mean that the 
accumulation would be less than 5% per year. Consequently, the value 
of "p" in the expression for g2 could decrease from 0.05 with time. 
On the other hand, current nominal traffic models would lead to 
between a 5% and 10% per year increase in the amount of US mass to 
orbit and some US and world traffic projections would give rise to 
increases in the accumulation of larger objects in orbit as high as 
20% per year. While such large increases do not seem historically 
justified, an upper limit of 10% increase per year, or p - 0.1, is not 
unrealistic. Any larger increases in the use of low Earth orbit would 
likely include different operational techniques which would invalidate 
assumptions used to express the design environment. 


Predicting the population not tracked by US Space Command is even more 
uncertain since we do not even have historical data to extrapolate. 
However, there are some indicators. Historically, the satellite 
fragmentation rate has increased with time, indicating that values for 
gl would increase with time faster than values for g 2 . However, 
actions are currently underway which should reduce the future 
satellite explosion rate. On the other hand, mathematical models 
predict that within the very near future, random collisions could 
become an important cause of satellite fragmentations. Under these 
conditions, the small debris population would increase at 
approximately twice the percentage rate of the large population, until 
a "critical density" of large objects is reached. This critical 
density corresponds to a value of g 2 between 10 and 100 (i.e., the 
tracked population is 10 to 100 times its 1985 total number) . At this 
time, values for would increase very rapidly with time, independent 
of values for g2 . 

The design environment assumes that the value of g^ increases at twice 
the percentage rate of g 2 . This could be expected if the satellite 
explosion rate continues to increase over the next decade or two. 
After this time random collisions would cause the rate to continue, 
independent of actions to reduce the explosion frequency. For values 
of p greater than 0.1, random collisions would become important in 
less than a decade, again consistent with the environment assumption. 
However, if the explosion rate is immediately reduced, and the current 
rate at which mass is placed into orbit does not significantly 
increase, then the design environment will predict fluxes for debris 
sizes smaller than 10 cm over the next 10 to 20 years which are too 
high by a factor of 2 to 10. 

III. Average mass density 

The average mass density for debris objects 1 cm in diameter and 
smaller is 2.8 grams/cm 3 . The average mass density for debris larger 
than 1 cm is based on observed breakups, area to mass calculations 
derived from observed atmospheric drag, ground fragmentation tests, 
and known intact satellite characteristics. 

This density has been found to fit the following relationship: 

p - 2 . 8 -d"° • 74 (4) 

IV. Velocity and Direction Distribution 

Averaged over all altitudes the non-normalized collision velocity 
distribution, i.e. the number of impacts with velocities between v and 
v + dv, relative to a spacecraft with orbital inclination i is given 
by the following equations: 


/(v) = ( 2 • v • v Q - V 2 ) • (G-e- ((v - a-v 0 )/(b-v „)) 2 + 
p. e -” v ' D -v 0 )/(E- Vo » 2 ) +H-C-(4-v-v 0 - v 2 ) 


( 5 ) 



where v is the collision velocity in km/s , A is a constant, and B, C, 
D, E, F, G, H, and v 0 are functions of the orbital inclination of the 
spacecraft. The values for these constants and parameters are as 
follows : 


A = 2.5 


B 


. 5 

- 0.5-0.01- (i-60) 

[ 0. 3 


i<60 

60<i<80 

i>80 


C 


'0.0125 i<100 

0.0125+0.00125- (i-100) i>100 


D - 1.3-0.01- (i-30) 


E- 0.55+0.005- (i-30) 


F 


0 . 3+0 . 0008 • (i-50) 2 
- 0.3-0.01- (i-50) 

0.0 


i<50 

50<i<80 

i>80 


G 


'18.7 

- 18.7+0.0289- (i-60) 3 
250.0 


i<60 

60<i<80 

i>80 


H - 1.0-0.0000757- (i-60) 2 


'7.25+0.015- (i-30) 
7.7 


i<60 

i>60 


When /(v) is less than zero, the function is to be reset equal to 
zero. An example for i - 30° is given in figure 4. 


The user may find it convenient to numerically normalize /(v) so that 


/'(v) 


/(v) 

oo 

f /(v) dv 
o 


( 6 ) 


When normalized in this manner, /' (v) over any 1 km/s velocity 
interval becomes the fraction of debris impacts within a 1 km/s 
incremental velocity band. Any average velocity moment may be defined 
as 


(7) 


oo 

v n — J v n */' (v) dv 
o 

The direction of impact can be approximated by using this velocity 
distribution and assuming that it results from the intersection of the 
spacecraft velocity vector and another circular orbit. That is, all 
velocity vectors will be in a plane tangent to the earth's surface, 
and will appear to be from a direction relative to the spacecraft 
velocity vector. The direction of the velocity vector is given by the 
relationship : 


cos $ = - ( 8 ) 

where 0 is the angle between the impact velocity vector and the 
spacecraft velocity vector, and v is the impact velocity. Since a 
spacecraft velocity of 7.7 km/s was used to calculate relative 
velocity, this velocity was used to determine the value of 15.4 
(2x7.7) given in equation 8 . 

A value for k (defined in section I) is found by integrating over the 
values of 0 that an oriented surface may be impacted. An example for 
i - 30° is given in figure 5, where the surface normal vector is 
located in a plane parallel to the Earth's surface, and has an angle 7 
to the spacecraft velocity vector. 


V. Uncertainty in Velocity and Direction Distribution 

The impact velocity and direction distributions are fundamentally 
functions of the orbital debris inclination distribution. The 
inclination distribution changes with time and altitude, and can 
change significantly as the result of a breakup at any particular 
altitude. Since the orbits of future breakups cannot be predicted, 
variables such as the altitude of the spacecraft are of secondary 
importance. Therefore, the most important variable is the inclination 
of the spacecraft. However, the velocity distribution will change 
with time and position in space. These changes could affect the 
average velocity from the distribution by several km/s. 

The fact that orbital debris objects are not in exactly circular 
orbits will introduce a small uncertainty for most velocities. As a 
result of the currently small eccentricities of these orbits, the 
actual direction of impacts are within 1 ° for most velocities derived 
from section IV. For low velocities (less than 2 km/s), the 
uncertainty in direction is much larger, with a significant fraction 
being more that 20° from the direction derived from section IV. This 
error in direction can be in the local horizontal plane or can appear 
as direction errors above or below this plane. High velocity impacts 
will almost always occur very near to the local horizontal plane and 
from the forward (down- range) direction; low speed impacts can occur 
from almost any angle ( 0 ° < angle < 180°) in the local horizontal 


plane as well as at considerable angles (0° < angle < 90°) out of that 
plane . 


VI. Flux Resulting from Possible Future Inadvertent Breakups 

The flux arising from the intentional or inadvertent fragmentation of 
an artificial earth satellite in low earth orbit (LEO) presents a 
hazard to other satellites. In the region of the breakup, an enhanced 
flux may be apparent for a considerable period of time, depending upon 
the altitude of the breakup, and the size and velocity distribution of 
the debris . 

The flux for a particle of mass m may be represented by the equation: 

F b - 1x10 “ 9 • <£ b • f * (M/m) (9) 

where F b is the flux of impacting fragments per square meter of 
surface per year, M is the total mass of the parent satellite, m is 
the mass of individual fragments in the same units as M, and f is the 
fraction of the total mass going into a fragment size characterized by 
m. This fraction may be derived from any differential number/mass 
distribution. The dimensionless quantity <£ b is a function of distance 
from the breakup altitude and the velocity of the ejecta from the 
center of mass; values for <£ b are given in Figure 6. 

To obtain values for 4> b , it was assumed that the breakup fragments 
were ejected in all directions from the center of mass of the parent 
object with a distribution of velocities. This distribution was 
assumed to have a "peak" or "most probable" velocity given by Vp, with 
the distribution linearly reducing to zero at 0.1-vh and 1.3-vh (i.e., 
on a number vs. velocity plot, the distribution is shaped like a 
triangle with the peak of the triangle at v b and a base range of 
0.1-v b to 1.3*v b ). Using this distribution of velocities, new orbits 
were calculated to obtain flux as a function of altitude. This flux 
distribution was then normalized and is depicted in Figure 6. 

The ejection velocity should not be confused with the collision 
velocity. The only time these two velocities would be identical is 
for the first few days following a breakup, and the object which 
fragmented is in the same orbit as the satellite at risk. However, 
the nodal crossing point of all orbits will precess at different 
rates, so that the collision velocity will increase with time. After 
a few years, the collision velocity would be close to the general case 
which depends on the orbital inclination. Inclinations greater than 
30° will yield collision velocities of 7 kra/s or greater. In general, 
the collision velocity will be similar to those given in section IV 
for most cases. 

The time for the flux to decay to e" 1 its initial value, or its "half 
life" H, for a 1 cm aluminum sphere and solar activity of S - 110, is 
given as a function of altitude in figure 7. When the breakup 
altitude is above the operational altitude, use the operational 
altitude to determine the half life. If the breakup altitude is below 
the operational altitude, use the breakup altitude to determine the 


half life. The half life is proportional to the particle mass- to-area 
ratio, so that the half life of other sizes can be derived. The total 
number of impacts resulting from a breakup is then 

N b - F b * A-H (10) 

where A is the surface area of a randomly oriented surface. Given the 
inclination of the breakup, both velocity and direction could be 
derived. 


VII. Discussion: An Example of a Future Breakup 

When a satellite breaks up in space, its size and velocity 
distribution are a sensitive function of the type of breakup. If it 
were a low intensity explosion, nearly all of the fragment mass would 
be in sizes larger than about 10 cm, and the most probable ejection 
velocity would likely be around 50 m/s. The fragments from a 
hypervelocity collision would include a significant fraction of mass 
with sizes less than 10 cm. However, the most probable velocities of 
these fragments would increase with decreasing size. Most of the 
fragments from a high intensity explosion could go into almost any 
preferred size, depending on the nature of the explosion. 

As an example, assume that half of the mass from a 1000 kg satellite 
goes into 1 cm fragments. Also assume that the satellite fragmented 
at an altitude of 600 km, and that the probable ejection velocity was 
150 m/s . The resulting flux of 1 cm fragments at 500 km would be 
5x10" 5 impacts/m 2 -yr. This is larger (by several factors) than the 
flux predicted at 500 km for 1995, given in section I. However, 
assuming no additional breakups occur, this larger flux will 
effectively last for only 3 years, as shown in figure 7. 
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Figure 1. The inclination scaling factor ^ (i) 
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Figure 2. A comparison of the predicted debris flux and the NASA 
SP 8013 meteoroid flux model (including Earth shielding) at 500km. 
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Figure 5. The ratio of flux on a spacecraft surface 
whose normal vector is oriented 7 degrees to the 
spacecraft velocity vector to the flux on a 
randomly-oriented surface, for an inclination at 30° 











